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SUMMARY 


Thirteen stainless- steel ring- stiffened cylinders were subjected to a pure 
bending load and heated rapidly until buckling occurred. For most of the cylin- 
ders the heating was not uniform around the circumference so that appreciable 
axial thermal stresses were present. Elementary thermal stress theory was found 
to be inadequate for the prediction of these thermal stresses, but a method was 
developed that would give satisfactory thermal stress results. By properly 
accounting for the thermal stress, the buckling load could be correlated with a 
theory for the buckling of an axially compressed ring- stiffened cylinder that is 
uniformly heated. 


INTRODUCTION 


A ring- stiffened cylinder used as a structural component in a missile or 
launch vehicle may undergo aerodynamic heating which, in general, is not uniform 
around the circumference or along the length. The resulting temperature varia- 
tion causes axial thermal stresses, and the presence of cooler rings causes cir- 
cumferential thermal stresses in the vicinity of the rings. These thermal 
stresses can cause reductions in the load carrying ability of the cylinder. The 
reduction of the axial compressive buckling load due to circumferential thermal 
stresses resulting from uniform heating is treated in reference 1, where a com- 
parison of theoretical results and experimental data is also given. If the cyl- 
inder is heated nonuniformly around the circumference, varying axial thermal 
stresses as well as circumferential thermal stresses are present. References 2 
and 3 present a theoretical analysis of the buckling behavior of a cylinder under 
a varying axial stress distribution. The results indicate that, theoretically, 
the maximum axial compressive stress at buckling is essentially the classical 
buckling stress for uniform compression. 

In order to determine the effect on buckling of combinations of axial and 
circumferential thermal stresses and load- induced stresses, an experimental 
investigation was conducted and the results are reported in this paper. Ring- 
stiffened cylinders were loaded by a pure bending moment and then heated non- 
uniformly until buckling occurred. An analysis of the axial thermal stress pres- 
ent in the cylinder is given in the appendix. The results of this analysis in 


combination with the theory of reference 1 is shown to give satisfactory correla- 
tion with the load- induced stress at cylinder buckling. 


SYMBOLS 


A = Et'R 6 

E'lLJ 

a. , b . 
in' in 

B = Et'R 2 
GtL 2 

c in> d in 


Fourier coefficients in expansion of shear flow in bay i 


Fourier coefficients defined by equations (Ak) 


E 
E * 

G 

I 

i 

L 

m 

n 

q. 

R 

T 

^max 


Young’s modulus of sheet material 
Young’s modulus of ring material 
shear modulus of sheet material 
moment of inertia of ring 
index for bay or ring 
length of cylinder between rings 

number of last bay in cylinder where first bay is zero (mil bays 
in a cylinder) 

general number of Fourier coefficient 

shear flow 

radius of cylinder 

temperature 

maximum cylinder temperature 


T 

AT 

AT 


temperature expressed as a function of 0 only and obtained by 
averaging temperature along bay length for a given value of 0 

change in T from initial or room- temperature value 

change in T from initial or room- temperature value 
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AT' 

t 

t’ 

u 

u in> u ' in 


average temperature rise of ring 
skin thickness 

skin thickness that is effective in carrying axial loads 
axial displacement 

Fourier coefficients in expansion of u^ 


x 

a 


axial coordinate 

coefficient of linear thermal expansion 


1 



constants defined by equations (A12), (A17), and (A20), respectively 


general notation for axial stress 

maximum axial compressive stress due to applied bending moment 


CJip 


a in> a 'in 


thermal stress 

Fourier coefficients in expansion of cx^ 


0 


circumferential coordinate 


TEST SPECIMENS 


The test specimens consisted of thirteen 19- inch- diameter cylinders having a 
nominal wall thickness of 0.030 inch with a resultant value of R/t of approxi- 
mately 300. All cylinders were of the same overall size with a total length of 


4^2 inches but there were two different ring spacings. Ten specimens had nine 

rings and three specimens had five rings with a resultant value of L/R of l/2 
and 1, respectively. The cylinders were fabricated by spotwelding. Because of 
limited sheet width it was necessary to have three longitudinal splices. Fabri- 
cation details of the specimens are shown in figure 1, and the ring spacing and 
cylinder wall thickness for each specimen are given in table I. 


The wall material was extra hard type 301 stainless steel and the Z- section 
rings were spun from type 304 stainless steel. The wall material was chosen so 
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that the variation of Young’s modulus with temperature would he small and the 
expected stress levels would he well within the elastic range of the material. 
The results then would not he a function of wide changes in material properties. 


TEST PROCEDURE 


Rapid Heating Tests 

A general view of the test setup and a detailed view of the heater installa- 
tion are shown in figure 2. One end of the cylinder is supported hy a heavy hack- 
stop; a pure bending load is applied at the other end hy a weight cage acting 
through a loading frame with pin- connected linkages. Rollers between the frame 
and support allowed lateral and longitudinal movement which minimized any extra- 
neous loads. The tip fixture was counterbalanced to eliminate any shear force on 
the cylinder. 

The cylinders were heated rapidly hy a 25- inch- long quartz lamp heater. For 
most of the tests the heater extended over approximately one- third of the circum- 
ference and was symmetrically located about the bottom extreme compression fiber 
of the cylinders, as shown in figure 2(b). In a few tests the heater extended 
around the entire cylinder in order to produce uniform heating around the 
circumference . 

The procedure for each rapid heating test was to apply a bending moment less 
than the room- temperature bending strength and then heat the cylinder at a rate 
of approximately 20° F/sec until buckling occurred. During each test temperatures 
at several ring and skin locations were recorded, and in some tests strains were 
measured with resistance- type strain gages. The vertical deflection of the tip 
fixture was continuously recorded in order to identify the instant of buckling. 

By varying the applied load in different tests, a buckling interaction curve for 
load- induced stress and temperature was obtained. 

Conventional room- temperature wire or foil gages were successfully used in 
regions where temperature did not exceed 175° F. The gage response due to tem- 
perature alone was not large compared with the overall strains being measured; 
therefore, the strains determined from these gages by using an appropriate 
temperature- response calibration curve were believed to be reasonably accurate. 

An attempt was made to measure strains in the heated portions of the cylinder 
with high- temperature strain gages. However, the gage response due to tempera- 
ture was large compared with the strains being measured in the temperature range 
of interest. The data obtained from the high- temperature strain gages were incon- 
sistent and generally could not be used as an accurate measure of strain in the 
cylinder. 


Room- Temperature Tests 

A limited number of tests conducted at room temperature were made in order 
to get an idea of the scatter to be expected in the data and to determine the 



reduction from the classical buckling stress for cylinders of the material and 
construction that were used. The room- temperature bending strength was determined 

for two cylinders having - = i and one cylinder having ^ = 1. Test procedure 

R 2 R 

was similar to that indicated for the rapid heating tests except load was applied 
in the absence of heating by means of a hydraulic jack instead of the weight cage* 


RESULTS AND DISCUSSION 


Buckling 

The results of the investigation are given in table I and are plotted in 
figure 3. For each specimen the maximum load- induced compressive stress is 
plotted as a function of the maximum cylinder temperature at time of buckling. 

The maximum cylinder temperature was taken as the temperature at the bottom or 
extreme compressive fiber of the cylinder and midway between rings in the bay 
nearest to the middle of the cylinder. This temperature was the actual maximum 
or nearly the maximum for all cylinders tested. 

It is of interest to compare the results at room temperature with predictions 
based on experimental cylinder buckling coefficients which are available from 
tests on aluminum cylinders (ref. k ) . It can be seen in figure 3 that the exper- 
imental room- temperature buckling data for the steel cylinders indicate a buckling 
coefficient less than the average but greater than the lower limit of the data 
on aluminum cylinders. There are not enough test data to establish whether any 
real difference in cylinder buckling coefficient exists for the two materials; 
however , any difference is probably small. 

If a cylinder were heated slowly and uniformly around the circumference so 
as to produce no thermal stress, the cylinder buckling stress would be reduced 
in proportion to the reduction in E. The amount of this reduction is indicated 
by the short- dash curves in figure 3 • Rapid but uniform heating causes circum- 
ferential thermal stresses in the vicinity of the cooler rings. The interaction 
of axial and circumferential stress results in a reduced buckling stress as shown 
by the circle test points in figure 3- The solid- line curve is based on the 
theory of reference 1 and tends to overestimate the test results somewhat. As in 
reference 1, it is necessary to modify the theoretical predictions based on the 
reduction in room- temperature cylinder buckling strength from the classical value. 
Inasmuch as this reduction had not been established for the type of cylinders 
tested, the results of the three room- temperature tests that were made were aver- 
aged and this value was used in determining the reduction in buckling stress from 
the classical value. The results of reference 1 are in the form of an interaction 
curve of axial compressive stress and the difference of the average skin and ring 
temperatures. In order to express the results in terms of it is necessary 

to know the variation of these two average temperatures with T max . This varia- 
tion is given in figure ^ where the average temperature of the interior bays and 
rings and also the difference in these temperatures are plotted as a function of 

T • The variation of a and E with temperature, which is also necessary for 
max 

the calculations, is given in figure 5. 
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For nonuniform heating, the load- induced stress at buckling is further 
reduced as evidenced by the square test points in figure 3. For this condition 
axial thermal stresses are present which might be expected to add directly to 
the load- induced stress so that buckling would occur when the total compressive 
stress (load- induced stress plus axial thermal stress) is equal to the buckling 
stress for uniform heating, as calculated from reference 1. Support for such an 
assumption is given in references 2 and 3, where buckling of a cylinder under 
varying axial compressive stress is shown to be relatively insensitive to the 
stress distribution but dependent on the maximum value of the compressive stress. 
When the load- induced stress at buckling is estimated by subtracting the axial 
thermal stress from the buckling stress for uniform rapid heating as calculated 
from reference 1, the long- short- dash curves in figure 3 result. Reasonable 
agreement is obtained except at the highest stress level where the results are 
somewhat unconservative for both the uniform and nonuniform heating tests. In 
the test of cylinder 6, failure occurred almost immediately after the heating 
started and the buckling temperature was much lower than was expected on the 
basis of results from the other tests. The values of the load-induced stress and 
buckling temperature were rechecked and no discrepancies were found. Therefore, 
a similar test was made on cylinder 7; this result was in agreement with the 
trend established by the other data. 

Calculation of thermal stress in cylinders involves more than the elementary 
thermal- stress computation, in which plane sections are assumed to remain plane, 
even if the internal region of the cylinder has a temperature distribution inde- 
pendent of length over many bays. The common practice, which applies to solid- 
section beams, of assuming that the effect of ends can be neglected at distances 
from the ends of the order of the beam depth is not a good assumption for cylin- 
drical shells. A method of calculating cylinder thermal stresses for a general 
temperature distribution is given in the appendix. This method was used in deter- 
mining the long- short- dash curve in figure 3- A detailed discussion of the ther- 
mal stress present in the test cylinders is given in the next section. 

The behavior of the specimens at buckling was typical of cylinders loaded 
in pure bending. Diamond- shaped buckles which extended to the vicinity of the 
neutral axis snapped in suddenly at the buckling temperature. Figure 6 is a 
photograph of cylinder 9 after buckling had occurred. 


Temperature and Thermal Stress Distribution 

A typical experimental temperature distribution along the bottom or most 
severely heated portion of the cylinder is indicated in figure 7 for a cylinder 

with i = i. The curves shown were faired through a number of data points. Tem- 
peratures were essentially symmetric about the center line and this symmetry was 
used in all calculations. Note the heat- sink effect of the rings and that the 
four interior bays received approximately equal heating. For the cylinders with 

^ = 1 the temperature distribution was essentially that which would be obtained 

I\ 

from figure 7 by removing the heat- sink effect at the even-numbered ring stations. 
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In the circumferential direction, temperature variations as shown in fig- 
ure 8 occurred midway between rings for cylinders heated on one side only. The 
curves shown are an analytical approximation of the data which was used in some 
of the thermal- stress calculations. 

The axial thermal stresses in the cylinders tested were calculated by using 
the analysis given in the appendix. A comparison between theoretical and exper- 
imental thermal stresses is shown in figure 9 where both longitudinal and circum- 
ferential thermal stress distributions are shown. The thermal stress is of most 
interest at the bottom (0=0) of the cylinder where the temperature is the 
greatest and where the load- induced compressive stress is a maximum. The bottom 
of the cylinder is the region assumed critical for buckling even though, as shown 
in the lower portion of figure 9(b), it is not the region of maximum compressive 
thermal stress. The compressive thermal stress is somewhat greater 20° away from 
the bottom because of the essentially constant temperature over a considerable 
portion of the cylinder circumference. However, the sum of the load- induced 
stress and thermal stress is a maximum at the bottom of the cylinder for most of 
the tests. At the higher temperatures the stress at the bottom of the cylinder 
is slightly less than the maximum stress. Nevertheless, the bottom is still 
probably the most critical region for buckling because the circumferential ther- 
mal stress which is also involved in the buckling is a maximum at the bottom of 
the cylinder and decays to zero in the unheated portion of the cylinder. 

As mentioned previously, accurate values of experimental thermal stress were 
not obtained over portions of the cylinder heated directly by the lamps. However, 
the comparison of measured and calculated stresses for the cooler areas does give 
an indication of the expected accuracy in calculating the thermal stress in the 
heated areas which are critical for buckling. In general, agreement between 
experiment and theory is good except that results of the calculations are greater 
than experiment in certain areas. The dashed curves in figure 9( a ) correspond to 
the elementary thermal stress in each bay and are seen to be considerably differ- 
ent from either the experimental or theoretical thermal stresses. 

One factor which was not accounted for in the theory was the effect of 
splices which are located at 0 = ±67tj and 0 = l80°. Along the splice at 

0 = l 80 ° which was the top of the cylinder and not in the region critical for 
buckling, the magnitude of the experimental stresses was, for the most part, less 
than that of the calculated stresses. Another factor which was not accounted for 
was the variation in the quantity aE. Above 600° F, aE decreases with tempera- 
ture so that the calculated thermal stresses are probably too high in this tem- 
perature range. The comparison of calculated and experimental buckling stresses 
(fig. 3(a)) also indicates that the theoretical thermal stresses are greater than 
the actual thermal stresses at the higher temperatures. 


CONCLUDING REMARKS 


The results of bending tests of rapidly heated, ring- stiffened stainless- 
steel cylinders have been presented. A few tests were made at room temperature 
and the experimental buckling coefficients were found to be essentially the same 
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or only slightly less than those which have been obtained for aluminum cylinders. 
The test results for uniform heating around the circumference could be correlated 
satisfactorily with a theory for the interaction of axial compressive stresses 
and circumferential thermal stresses. When the heating was nonuniform, substan- 
tial axial thermal stresses were present. These stresses, which were consider- 
ably different from those predicted by elementary theory, showed reasonable agree- 
ment with a theory that takes into account radial deflections of the skin and 
rings. The load- induced stress at buckling for the nonuniform heating tests could 
be satisfactorily predicted by subtracting the calculated axial thermal stress 
from the predicted buckling stress for uniform rapid heating. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August 29, 1962 . 
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APPENDIX 


AXIAL THERMAL STRESS IN A RING- STIFFENED CYLINDER 
WITH AN ARBITRARY TEMPERATURE DISTRIBUTION 


The elementary theory of thermal stress can be used satisfactorily for 
solid-section beams if the interior portion has a temperature distribution essen- 
tially independent of length. In the vicinity of a free end the thermal stress 
must decay to zero. However, this condition usually is of little concern because 
only the maxi mum thermal stress which occurs away from the ends and which can be 
calculated by elementary theory is generally required. The elementary theory, 
when applied to a ring-stiffened cylinder, does not account for radial deflec- 
tions of the rings or cylinder wall. Inasmuch as the dimensions of most ring- 
stiffened cylinders are such that these deflections cannot be ignored without 
introducing appreciable inaccuracies in analytical results, the elementary theory 
must be modified so as to satisfy continuity between rings and sheet as well as 
the equations of statics. The following development yields self-equilibrating 
stress distributions which can be superimposed on the elementary thermal stress 
distribution so as to maintain the continuity of the cylinder. The method is 
essentially the same as that of reference 5 where self-equilibrating stress 
distributions were found that could be superimposed on the stresses obtained 
from the elementary theory of torsion and bending to yield the correct stress 
distribution for a cylinder subjected to lateral loads. 


General Theory 

For purposes of analysis the cylinder is divided into bays small enough 
that the temperature in each bay is a function of 0 only or 

Tq = T ± (0) (Al) 

where i denotes the bay number of the cylinder which has a total of m + 1 
bays^ The bay notation and coordinate system are shown in figure 10. The bar 
in T^(0) denotes that at a given value of 0 an average temperature over the 

bay length is used. The most satisfactory arrangement is obtained when the ring 
spacing is small enough that equation (Al) is a reasonable approximation of the 
temperature for bay widths equal to the ring spacing. If the temperature varia- 
tion in the axial direction is too large between rings, several bays may be 
required between rings. In this case the bays between the actual rings are 
assumed to be separated by a fictitious ring having a moment of inertia based on 
the effective sheet material in each bay. With this assumption it is seen that 
the bay length and ring spacing are not necessarily equal. 

The thermal stress in the cylinder may be obtained by a process similar to 
that used in elementary thermal stress problems. If each bay is rigidly con- 
strained, the resulting stress distribution is given by 
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a ± = -a^ AT.(0) (A 2) 

In each bay this stress system is independent of x^, but for equilibrium the 
constraints at the end of each bay are loaded by forces with intensity a^t \ 
where a ^ is given by equation (A2). If a loading intensity of -cr^t'^ is 

applied to the cylinder at the ends of each bay, the forces on the constraints 
are removed. Therefore, the constraints are no longer needed to maintain the 
cylinder in equilibrium. Thus, the thermal stress is calculated by adding the 
stress given in equation (A2) to the stress system obtained from loading the 
cylinder at the end of each bay with stresses equal but of opposite sign to those 
given in equation (A2). 

Solution of stress-distribution problem .- In solving this stress-distribution 
problem it is convenient to express the applied stresses at the end of each bay 
as a Fourier series 


o 1 E i AT ± = ^ (c in cos n0 + d ±n sin n0 ) (A3) 

n=0 

The Fourier coefficients c^ n and d^ n may be determined as 



p 2it 

d in = “ / sin n0 d0 

* J 0 J 


The stress distributions corresponding to c_^, c i3/ ^ii axe a un ^^ orm 

axial stress distribution and two bending-moment distributions about two perpen- 
dicular axes. These stresses do not vary with Xj_ and are seen to yield the 

elementary thermal stress solution for bay i when added to equation (A 2). The 
remaining terms (n > l) represent a self-equilibrating stress distribution which 
decays with increasing distance away from bay i. 

The solution of the stress-distribution problem for n > 1 can be obtained 
by an extension of the method presented in reference 5 . In the theory of refer- 
ence 5 it is assumed that the strain in the circumferential direction is zero. 
All material which resists axial stresses (including longitudinal stringer 
material) is represented by an effective thickness t 1 ^, whereas the actual 

thickness t is assumed to resist only shear stresses. By using these 
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assumptions, the self -equilibrating shear flow and direct stress in bay i and 
the associated axial displacements can be written as follows (see ref, 5 )' 


00 

q ± (0) - sin nS?S + bin cos n ^) 

n=2 

OO — 

Oi( x i>0) = ^ a in( x i) cos n S i + ff 'in( x i) sin n 0 
n=2 *- 

OO — 

u^x^jZl) = £ u in (x ± ) cos n0 + u* in (x i )sin n0 
n=2 L 


where 


,(o) 


/ \ 

a in( x i) = - ^ a in + 0 in 


0 ’in( x i) = S 7 b in + 0 ’in(°) 




and 




u in( x i) = - 


u 'in ( x i) ~ 


nx 




2E i Rt i 


aj ^ + 


E, 


x i + u in<°) 




mcj 


2E i Rt ' ± 


b in + 


°in(°) 

E. 


x i + “'in* 0 ) 


J 


(A5) 

(a6) 


(A7) 


(AS) 


(A9) 


The complete stress distribution can be obtained once the Fourier coefficients 
a^ n and b^ n sure determined. The solution for a^ n and b-^ n for the problem 

considered herein is the same as that in reference 5 except that the equations 
for equilibrium of forces at a ring (see fig. 10) are changed as follows; 

[ 0 i-l,n( L i-l) " c i-l, n] b 'i-l = [ 0 in (O) ‘ 

[ c 'i-l,n( L i-l) - d i-l,n] b 'i-l = [ 0 'in (o) “ 

In reference 5 the thicknesses t\ are omitted in error and there are no 
applied in-plane stresses so that c^ n and d.j_ n do not appear. 
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From equation (A3) it can be seen that the coefficients ci n correspond 
to symmetrical loadings and the coefficients dp n correspond to antisymmetrical 
loadings. Since the solution for either loading is essentially the same, only 
the solution for the symmetrical terms, corresponding to the heating condition 
for the cylinders reported herein, are considered in the following development. 

A recurrence formula for the symmetrical shear-flow coefficients a-j_ n is 

obtained as in reference 5 in the form 




E'i-l 1 !-!' 


1 + 


E 'i-I r i -1 , ^ B i -1 “ nC 


E'il. 


6 A i-l^ / E 'i I i 


+ 


1 + 


E’iIi E 'il-tl f i T 
ii i i\ hi i+1 i i+1 


E 'i I i , 6B i " 
6Aj7 


v 2 + v U n 


E f I 2E ’ a I n * A,* 7 

* i+l 1 i+l 1 1 ^ L/ 


a i+l,n 


E 'i + i i i+i e ' i+i^i+i 


+ fWia + 6 B i*i ' 


P I T 

i+ 2 x i +2 


6A i+i? 


+ a i+2,n(;r; 7 

\ E 1+2-4 +2 > 


2R 5 7 


t 2 

v l L i-l 


^♦1 


± l-x / . , \ y i + x / , , , v 

i 77 rin 1 in " c i-l,n t i-l,n) + 7 77 ( c i+l t i+1 " c in t i) 

i-l 1 i-l v ' E i+l t i+l v 1 


(All) 


where 


(2 ^ i < m - 2) 


. i / E i t 'i L i-i + hzl 
V 1 ' L 2 U_ n t',_ n L, L, 


i-1 


“i-l 1 ’ i-1 i H+1j 


Vo * 


V-, = 


i/ 2L i VA-I 


Lf\ L i+l 


%-!* 'i-l L i+l 


^’i 
E i+l t 'i+l/ 


E.t'.L. , 
l i i-l 


) (A12) 


5 


+ 1 


Vi, = 


1 / L i + E i t 'i L i-l 
L i \ i+l 


1 2 

l 1 i L i-l 


Wi-A+l 


E i-l E i+l t 'i-l 1 'i+l L i 


^'i > 

E i+l t 'i+l) 


J 


12 



and 




1 



E ± t ’jR 6 

% = 5 

E’iliLi 

Ejt \R 2 

B i = 2" 

GitiLi 




J 


(A13) 


B oundary conditions .- Equation (All) applies to any interior bay i when 
there are at least two adjacent bays on both sides of the given bay. This con- 
dition is not satisfied at the ends of the cylinder and special equations are 
necessary which are determined by the boundary conditions. In all cases one 
boundary condition is that the radial deflection of the cylinder wall is equal 
to the radial deflection of the ring at the end of the cylinder. When there is 
no axial stress on the end cross section, which is sometimes referred to as a 
free end, an additional boundary condition is 

W°> - <=0n ‘ 0 (M 


The following equations axe obtained as a result of these boundary conditions: 


and 




A 1 

/ ( E- 0 I 0 t 6B 0 - n A 

Ap A*n 2 

+ ^ ~ 

E 'o I o 

l E'xii 6 v y 

E, i I i 


a ln 


r ' 2 n + E 'i I i + 6B i ~ n ^ ) + A i 


E'jlA E '2 I 2 6A 1? / E 'l I l 


Me' 2 i 2 


p ( A l L O c On + ?v 2 L l c In | 


2r5 7 \ e o 


(A15) 


A^n' 


0n\ E 1 3_Ix 2E’ 0 XoAo7 


ft ln 


*2 L + iiii + 6B 1 - n !) + A_ 

E ' x 1 ! \ E '2 I 2 6A 1 7 / E, 2 I 2 




2E'i l i A 1 r 


2r5? 


E- 2 I 2 E-2l 2 


*2 + \ L + k 2 x 2 + on 2 


E 1 pi o 6B p — n 


e '5 I 5 


6^7 


Me' 3 i 3 


fin + \ L 2 


“■1 E- 


~r( c 2a t '2 - C ln t 'l) 
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( Al6) 
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where 


A 1 


1 

L 0 L 1 



*3 



. Wo\ 
Ei t \J 




1 

L 1 L 2 


1 




A = 1 / E O t ' o L 2 + E O t ' Q \ 

Lq \ E 2 t ' 2 L 1 E l tf l / 


(A17) 


^6 


3 'l L 2 

if \ E 2t '^la 



If the end of bay zero is clamped, the axial displacements due to an axial 
force and a bending moment may be allowed, but the displacements due to the self- 
equilibrating force system are constrained so that Uq(O, 0) is equal to zero 

which implies UQ n (o) is equal to zero. The following equations are obtained 
for a clamped end: 


and 


_!1_ + ^ L + iVo + 6B o - f \ + ^ ng 
E'l 1 ! E'qIoV E'ili 6Aq? J 2E' 0 I 0 Ao7 


4 1 
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where 
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L 1 
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^2 = 




E 0 t ' 0 
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The notation used in equations (All) to (A20) is essentially that of reference 5 
except for the presence of a different Young's modulus for each bay and ring. 
Boundary conditions for the opposite end of the cylinder, which is denoted by 
bay m, may be obtained from equations (A15) to (A20) by a suitable change of 
subscripts. The subscripts which refer to bay numbers 0, 1, and 2 are replaced 
with the subscripts m, m - 1, and m - 2, respectively. Subscripts 0, 1, 2, 
and 3 which refer to ring numbers and occur in the term E'^I-^ are replaced by 

m + 1, m, m - 1, and m - 2, respectively. It is also necessary to change the 
algebraic sign of the right-hand side of equations (AI 5 ), (Al 6 ), (Al 8 ), and 
(A19). The condition of zero radial deflection of the cylinder wall at the ends 
can be specified by allowing the appropriate ring moment of inertia (I 0 or I m+ j} 

to approach infinity. For the limiting case of rigid rings where all moments of 
inertia Ij^ approach infinity, the governing equations are greatly simplified. 


The application of equation (All) and the appropriate boundary conditions 
results in m + 1 equations for the determination of the m + 1 different 
values of a^. If the stress at the end of a cylinder is specified as zero, 

the axial stress distribution can be determined from the coefficients apn an( ^ 

15 



the use of equations (a 8 ) and (A10). If the axial displacement is specified as 
zero at the ends, the axial stress distribution is determined from the coeffi- 
cients a^ n and equations (A 8 ) to (A10). Solutions for unsymmetrical loadings 

can be obtained by replacing aq n with -b in in equations (All), (AI 5 ), (Al 6 ), 
(Al 8 ), and (A19). 

The axial thermal stress at any ring can now be written for the general 
case as 


drjq(O) = + c^ cos 0 + d^ sin 0 

00 

^ [5i n (°) cos n 0 + n 0J 


(A21) 


n=2 


There is a linear variation of ov^ between each ring. Note that the first 

four terras on the right-hand side of equation (A2l) represent the elementary 
thermal- stress solution for each bay. 


Application of Equations to Test Cylinders 

The equations presented thus far are very general in that each bay length 
and ring size may be different. Also material properties may be different for 
each bay and ring. However, the material properties are assumed to be constant 
around the circumference for a given ring or bay; this means some average or 
effective value of E must be assigned to each ring and bay if there is a 
circumferential, variation in E. The manner in which E appears in the recur- 
rence formula and boundary conditions is in the ratio of E for the material 
in one bay or ring to E for the material in a nearby bay or ring. These equa- 
tions involve, at the most, five adjacent cylinder bays; if E does not vary 
greatly within these bays, the ratio of the values of E for different bays may 
be assumed to be unity. With this assumption the actual (or effective) value of 
E is not required for the solution of equations (All), (A15), (Al 6 ), (Al 8 ), 
and (A 19 ). This assumption is made in the calculations for the test cylinder of 
this report. The calculation of c-j_ n and dj_ n also involves E but only in 

the product clE which tends to be more constant with temperature than either a 
or E. A constant value of 0.24 ksi/°F was assumed in all calculations. 

The axial thermal stress for the test cylinders may be determined by the 
following procedure. A schematic diagram of a typical test cylinder is shown 
in figure 11 which also includes the dimensions and ring moments of inertia 
necessary for the calculation. All cylinders were of the same overall size but 
there were two values of ring spacing. The idealization shown in figure 11 has 

equal bay lengths corresponding to i For = 1 the same bay length is 

n 2 n 

used but the even-numbered ring locations, which correspond to midway between 
rings in the actual cylinder, are assigned a moment of inertia equal to that of 
the wall material one bay wide. The end bays are actually shorter than the other 


16 



bays and terminate at a 5 / 16 - inch adapter ring which is fastened to a very heavy 
attachment ring. (See fig. 1.) The whole assembly is then bolted on one end to 
a very stiff tip fixture and on the other end to a backstop. Because of this 
stiff end attachment, the end bays were assumed to be clamped and to have a 
length equal to that of the other bays. The heating is assumed symmetrical about 
0=0; therefore, only the equations involving a-j_ n are needed. Because of 

symmetry about ring 5 > a*_ n * ~ a 4 n> a 6 n = ~ a 3 n' anc ^ c ^ n * c 4 n> thus, equations 

for only five bays have to be considered. With these assumptions and boundary 
condition, the equations determining aj_ n can be obtained from equations (All), 
(Al 8 ), and (A19) in a simplified form as follows: 


a 0n( 2J|A l? + l8B + 5n2 ) - a ln( 24 V + 6 V + 6B - n2 ) + a 2n( 6 V) = ^P( c ln ‘ c 0n) 

-a 0n( l8A l 7 + 6b " n2 ) + a i n ( l8A l 7 + i8a 2 y + 12B + ^ n2 ) 

- a 2n (6A 1 7 + I8A27 + 6b - n 2 ) + 


3Rt'n 

= ~ "On 


( c 2n _ c 0n) 

a 0n (6A 1 7) - a ln (6A 1 7 + 18a 2 7 + 6 b - n 2 ) + a 2n (l8A 1 7 + l8Ag7 + 12B + 4n 2 ) 


(i 8 a x 7 + 6Ag7 + 6B - n 2 ) + a J+n (6A 2 7) 


a 3n 


3Rt'n, x 

— ( C ?n - C ln) 


In 


[GP^y) - a 2n (l8A 1 7 + 6^7 + 6 b - n 2 ) + a 3n (l8A 1 7 + 18^7 + 12B + 4n 2 ) 
- a lfn (l2A 1 7 + 18A 2 7 + 6B - n 2 ) 

Kn ‘ C 2n) 


3Rt'n, 


a 2n( 8A l 7 ) " a 3n( 8A l 7 + 2l+A 2 7 + 8b " r2 ) + a 4n(^ A l 7 + 2AA 2 7 + l8B + 3 r2 ) 


3n 

^P( C 4n - c 3n) 


) (A22) 
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In equations (A22) the subscripts on L, B, and t' have been dropped inasmuch 
as there is no change in these quantities with bay number. For cylinders with 

— = 1 there are two independent values of Ap given by Ap and A 2 . 

R 


The calculation of the coefficients c in by use of equations (A4) requires 

some special consideration. The final determination of the thermal stress often 
involves small differences of large numbers. An error of less than 1 percent in 
determining cp n in some instances can easily cause errors greater than 10 per- 
cent in thermal stress. Thus, it is especially important to obtain an accurate 
consistent set of values for c in . Graphical or numerical integration of equa- 
tions ( Ak ), unless done very accurately, may lead to serious errors in the final 
result. For this reason, an analytical expression for the temperature distribu- 
tion was obtained that yielded an analytical expression for c in . Thus, c in 

necessary and all values would be consistent 
the temperature rise was applied to all bays 

"S 

0 + 0.05 sin 50 (0° ^ 0 £ 60 °) 

0 ( 60 ° ^ 0 £ 8o°)> (A 23 ) 


could be determined as accurately as 
with each other. The expression for 
and is given by 

ZTp 

— = cos 

ATi(O) 

ATp 

— = cos 

Zt.(o) 


Alp 


ATpCo) 


= COS 


80 ° e -°-°623(0-8o o ) (go 0 % 0 ^ 180 °) 


j 


This temperature distribution was taken to be symmetrical about 0-0. The bars 
denote an average temperature rise over a bay length. This average temperature 
rise in a given bay was essentially a constant tiroes the temperature rise at 
midbay. Thus, the expressions on the right-hand side of equations (A23) apply 
equally well to the ratio AT i ^AT i (0). The curves of figure 8 were calculated 

from equations (A23) by replacing ATp^AIp(O) with ATp0ATp(O), and the agree- 

ment between measured and calculated values is indicative of the accuracy of 
equations (A23) for both the average and midbay temperatures. Values of ATp(O) 

and Afi(0), which are necessary to calculate the temperature distribution, were 
obtained for each bay from the measured temperature distributions and are given 
in figure 12 as a function of the maximum cylinder temperature. 


After the coefficients a in have been determined, cr 0rl (0) can be obtained 
from equations (a 8) to (AlO) with the use of the conditions u 0n (0) = u^ n (0) = 0: 


W°> * J 

[2 nEt'( 9a 0n + 7a ln + 5 a 2n + >3n * a 4n) 



Gin + ° 2 n + c 3 n + c 4n ' Uc 0 nj] 

( A2k) 
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The remaining coefficients cr in (0) are determined from the first of equa- 
tions (A8). The thermal stress in the test cylinders was calculated by using 
equation (A2l); a comparison of the calculated and experimental thermal stresses 
is shown in figure 9* Terms up to n = 10 were used in the calculations. This 
number of terms resulted in suitable convergence except for the stresses in 
bay 2, which are estimated to deviate less than 5 percent from the value that 
would be obtained by taking many more terms. This behavior is believed to be 
caused by the steep longitudinal temperature gradient in this bay. The assump- 
tion of clamped ends is supported by the agreement of experiment and theory in 
the vicinity of the ends. A characteristic of cylinder behavior which is con- 
trary to that observed in solid-section beams is the slow decay of the thermal 
stress away from the heated area. (See fig. 9(a).) The stress distribution is 
affected by ring stiffness as illustrated in figure 13 where the axial thermal 
stress along the bottom of the cylinder is plotted for three values of ring stiff- 
ness. The ring stiffnesses used were: completely rigid = 00 ), the actual 

stiffness as measured by the calculated ring moment of inertia, and one-half the 
calculated stiffness. Little difference is seen between the curves for rings of 
finite stiffness but the stress distribution obtained by assuming rigid rings 
exhibits a much greater rate of stress decay. Even though the rings used were 
reasonably stiff by most standards, their behavior for thermal stress was far 
from that of a rigid ring. 
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TABLE I 

RESULTS FOR BENDING AND HEATING TESTS OF RING- STIFFENED CYLINDERS 


Cylinder 

in. 

l/r 

a x> 

ksi 

T 

- L max* 

Op 

Type of 
heating 

i 

0.0506 

1/2 

35-1 

Room 

None 

2 

.0301 

1/2 

38.1 

Room 

None 

3 

.0305 

1 

4 o .4 

Room 

None 

4 

.0305 

1/2 

29.4 

380 

Uniform 

5 

.0503 

1/2 

24.4 

800 

Uniform 

6 

.0306 

1/2 

29.4 

120 

Nonuniform 

7 

.0306 

1/2 

29.4 

250 

Nonuniform 

8 

.0307 

1/2 

24.1 

450 

Nonuniform 

9 

.0304 

1/2 

19.1 

625 

Nonuniform 

10 

.0305 

1/2 

16.9 

833 

Nonuniform 

li 

.0304 

1/2 

11.85 

910 

Nonuniform 

12 

.0301 

1 

25.1 

400 

Nonuniform 

13 

.0311 

1 

19.2 

635 

Nonuniform 
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(a) General view of test setup. L-6I-6326.I 



(b) Detail view of heater installations. L-6I-6327 
Figure 2.- Method of loading and heating test cylinders. 
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Figure 6.- Buckle pattern for cylinder 9- L-60-3^4 
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Eq uilibrium of axial forces at a ring 

Figure 10.- Cylinder coordinate system with bay and ring numbering system. 



I o~ 1 1 0 = 00 

For 

L/R = 1/2 


t'=t = 0. 0306 i n . 


Ij = 0.0295 inf 


R = 9 . 56 i n . 

For 

L/R = 1 




Ij = 0.0295 inf 

( i = odd ) 



I: = 0.00001 13 inf 

( i = even ) 


Figure 11.- Idealization of test cylinders for analysis. 
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(b) Average temperature rise. 

Figure 12.- Variation of temperature rise at bottom of cylinder with maximum 

cylinder temperature. 












